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Real complex systems are inherently time-varying. Thanks to new communication systems and 
novel technologies, it is today possible to produce and analyze social and biological networks with 
detailed information on the time of occurrence and duration of each link. However, standard graph 
metrics introduced so far in complex network theory are mainly suited for static graphs, i.e., graphs 
in which the links do not change over time, or graphs built from time- varying systems by aggregating 
all the links as if they were concurrent in time. In this paper, we extend the notion of connectedness, 
and the definitions of node and graph components, to the case of time-varying graphs, which are 
represented as time-ordered sequences of graphs defined over a fixed set of nodes. We show that the 
problem of finding strongly connected components in a time- varying graph can be mapped into the 
problem of discovering the maximal-cliques in an opportunely constructed static graph, which we 
name the affine graph. It is therefore an NP-complete problem. As a practical example, we have 
performed a temporal component analysis of time- varying graphs constructed from three data sets 
of human interactions. The results show that taking time into account in the definition of graph 
components allows to capture important features of real systems. In particular, we observe a large 
variability in the size of node temporal in- and out-components. This is due to intrinsic fluctuations 
in the activity patterns of individuals, which cannot be detected by static graph analysis. 

PACS numbers: 89.75.Hc,89.75.-k,89.75.Fb 
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Time-varying graphs are a natural model for 
networked systems in which the relationships 
among nodes are intrinsically dynamic and fluctu- 
ate over time, where links appear and disappear 
at specific points in time and are often recurrent. 
Here we extend the concept of connectedness and 
the definitions of node and graph components to 
the case of time-varying graphs, we prove that 
finding strongly connected components in time- 
varying graphs is an NP-complete problem and 
we also report the results of component analy- 
sis performed on three real time- varying systems. 
This analysis confirms that the classical aggregate 
representations of networks evolving over time 
wash out most of the richness of the original sys- 
tems. In particular, static graph erroneously flat- 
ten down fiuctiations in the size of in- and out- 



componens of nodes, and tend to substantially 
overestimate the actual size of the connected com- 
ponents of the graph. 



I. INTRODUCTION 

Complex network theory has proved to be a versatile 
framework to represent and analyze biological, social and 
man-made complex systems [l|, ■ Typically, a complex 
system is inherently dynamic. Social interactions and hu- 
man activities are intermittent [H-lfl, the neighborhood 
of individuals moving over a geographic space evolves 
over time [7j9|, links appear and disappear in the World 
Wide Web 101 ] . in patterns of interactions among genes 
from microarra y ex periments fill Il2l l and in functional 
brain networks ll!; 



. In all these networks, time plays 
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a central role: links exist only for certain time periods, 
and are often recurrent. Despite this fact, most of the 
classic studies in complex networks theory are based on 
the analysis of the topological properties of static graphs. 
These are graphs in which the links do not change over 
time, or graphs built from time- varying systems as the 
result of the aggregation of all interactions, as if these 
were all concurrent in time. The evolution of linking 
patterns over time, when considered, has been usually 
studied by creating a series of graphs, each graph con- 
taining all the links appeared in a certain time interval. 
Then, each standard graph metric has been evaluated 
for the static graph obtained at each time window, and 
plotted as a function of time [HI, [l6[ ■ Today, thanks to 
recent technological developments, for the first time we 
have the opportunity to study large social and biological 
networks with precise temporal information on the ap- 
pearance, duration and frequency of links among a set of 
nodes. Many other similar databases will be produced 
in the near future, at an ever increasing rate. These 
data sets demand for new network measures and mod- 
els that can take account of the richness introduced by 
detailed temporal information. Some recent works have 
analyzed large interconnected systems with fluctuating 
interactions [l7H19| , and some graph measures have been 
already extended to the case of graphs in which connec- 
tion patterns evolve over time [20[. More recently, some 
works have also studied the onset of synchronization in 
populations of agents interacting throu gh t ime-evolving 



topologies [2l|. In previous works |22H24| it has been 



shown that a static analysis of aggregated graphs is not 
able to capture the real dynamic behavior and time cor- 
relations of complex networks evolving over time. Since 
in this type of analysis the temporal ordering of links 
is ignored, the number of available links at each time is 
overestimated and, therefore, the actual lengths of walks 
and paths is underestimated. In particular, it has been 
found that edge causality plays an important role on the 
dynamics of many processes that occur on networks, such 
as disease spreading [25l - |27j . There have been several 
attempts to incorporate temporal information in the de- 
scription of complex networks [25|, [H, [29| , and different 
definitions of time- varying graphs have been proposed so 
far [32l-f34]]. At the same time, some basic concepts of 
complex network theory, such as temporal walk s f30L l3lj , 
path length and distance [3214341 ]. centrality [23j] have 
been extended to the case of temporal networks. 

In this paper, we focus our attention on two important 
concepts in graph theory, namely those of connectedness 
and connected components of a graph, and we generalize 
them to the case of time- varying graphs. These concepts 
have been thoroughly used to study the reachability of 
pairs of nodes in static complex networks (35j . and to 
characterize the resilience of networks to attacks j36j . 
Here we show that node connectedness and connected 
components play a central role in time-varying graphs, 
and can reveal interesting details on the real structure of 
the network, that usually remain hidden to an analysis 



based on static graphs. 

The paper is organized as follows. In Section [ill wc 
briefly review the concepts of connectedness and com- 
ponents in static graphs, while in Section IIIII we extend 
them to the case of time-varying graphs. In particular, 
we define the temporal in- and out-component of a given 
node, and we give the definition of weakly and strongly 
connected components of a graph. In Section IIVI wc 
show that the problem of finding components in a time- 
varying graph can be mapped into the maximal-clique 
problem for an opportunely constructed graph, which 
we call affine graph. An affine graph is a static graph 
which incorporates all the information on the temporal 
reachability of pairs of nodes, and is a useful tool to ana- 
lyze the components of the corresponding time-varying 
graph. Thanks to this mapping, we also prove that 
finding strongly connected components in time- varying 
graphs is a NP-complete problem. Finally, in Section [V] 
we present the results of temporal component analysis 
on time-varying graphs constructed from three different 
data sets of human interactions. 



II. COMPONENTS IN STATIC GRAPHS 

Let us consider a graph G with N nodes and K links. 
From now on we will refer to it as to a static graph. We 
will consider the case of undirected and directed static 
graphs separately. An undirected static graph G can be 
represented by a symmetric adjacency matrix, a N x N 
matrix A whose each entry a%j is equal to one if and 
only if there is a link between i and j, and is equal to 
zero otherwise. In order to define graph components, 
we need to introduce the concept of connectedness, first 
for pairs of nodes, and then for the whole graph. Two 
nodes i and j of an undirected graph G are said to be 
connected if there exists a path between i and j. G is said 
to be connected if all pairs of nodes in G are connected, 
otherwise it is said to be unconnected or disconnected. 
A connected component of G associated to node i is the 
maximal connected induced subgraph containing i, i.e., 
the subgraph induced by all nodes connected to node 
i. If an undirected graph is not connected, it is always 
possible to find a partition of the graph into a set of 
disjoint connected components. It is straightforward to 
prove that this partition is unique. 

A directed static graph G is described in general by 
a non-symmetric adjacency matrix, a N x N matrix A 
whose each entry aij is equal to one if and only if there is a 
directed link from i to j, and is equal to zero otherwise. 
Defining connectedness for pairs of nodes in a directed 
graph is more complex than in an undirected graph, be- 
cause a directed path may exist through the network from 
vertex i to vertex j, but this does not guarantee that any 
path from j to i does actually exist. Consequently, we 
have two different definitions of connectedness between 
two nodes, namely weak and strong connectedness. In 
particular, we can define the weakly and the strongly 



connected components of a directed graph as follows [37| . 
Two nodes i and j of a directed graph G are said strongly 
connected if there exists a path from i to j and a path 
from j to i. A directed graph G is said strongly con- 
nected if all pairs of nodes are strongly connected. 
A strongly connected component of G associated to node 
i is the maximal strongly connected induced subgraph 
containing node i, i.e., the subgraph induced by all nodes 
which arc strongly connected to node i. A weakly con- 
nected component of G is a component of its underlying 
undirected graph G w , which is obtained by removing all 
directions in the edges of G. Two nodes i and j of G 
are weakly connected if they are connected in G u , and a 
directed graph G is said to be weakly connected if the un- 
derlying undirected graph G u is connected. Hence, the 
components of a directed graph can be of two different 
types, namely weakly and strongly connected. It is also 
useful to review the definitions of components associated 
to a node of a directed graph. We have four different 
definitions: 

1. The out-component of node i, denoted as OUT(i), 
is the set of vertices j such that there exists a di- 
rected path from i to j, Vj. 

2. The in-component of a node i, denoted as IN(i), is 
the set of vertices j such that there exists a directed 
path from j to i,Wj. 

3. The weakly connected component of a node i, de- 
noted as WCC(i), is the set of vertices j such that 
there exists a path from i to j, Vj in the underlying 
undirected graph G u . 

4. The strongly connected component of a node i, de- 
noted as SCC(i), is the set of vertices j such that 
there exists a directed path from i to j and also a 
directed path from j to i,Vj. 

We have implicitly used the last two concepts for the def- 
initions of weakly and strongly connected components of 
a directed graph given above. In fact, the property of 
weakly and strongly connectedness between two nodes 
is reflexive, symmetric and transitive, i.e., in mathe- 
matical terms, it is an equivalence relation. Therefore, 
it is possible to define weakly and strongly connected 
components of a graph by means of the weakly and 
strongly connected components associated to its nodes: a 
strongly (weakly) connected component of a node is also 
a strongly (weakly) connected component of the whole 
graph. Conversely, the definitions of out-component and 
in-component of a node are not based on equivalence re- 
lations. In fact, the symmetry property does not yield: 
i G OUT(j) does not imply j € OUT(i). This means that 
out- and in-components can be associated only to nodes, 
and cannot be directly extended to the entire graph. In 
practice, we cannot partition a graph into a disjoint set 
of in- or out-components, while it is possible to identify 
a partition of a static graph into a disjoint set of weakly 
or strongly connected components. However, the in- and 
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FIG. 1. A directed graph can be partitioned into a set of 
disjoint weakly connected components (in yellow). Further- 
more, each of these components has a rich internal structure, 
as shown for the GWCC. 

out-components of the nodes of a graph can be used to 
define the strongly connected components of the graph. 
From the above definitions, we observe that i G OUT(j) 
if and only if j G IN(i). Furthermore, we notice that i 
and j are strongly connected if and only if j G OUT(i) 
and, at the same time, i G OUT(j). Or cquivalently, if 
and only if j G OUT(i) and j G IN(i). Therefore, the 
strongly connected component of node i is the intersec- 
tion of W(i) and OUT(i). 

We are now ready to describe the rich interplay among 
the various concepts of connectedness in a directed static 
graph in detail. In the most general shown in 

Fig. [T] a directed graph can be decomposed into a set 
of disjoint weakly connected components. In a large 
graph, one component will be larger than all the oth- 
ers. This component is usually called the giant weakly 
connected component GWCC of the graph. If we treat 
each link in the GWCC as bidirectional, then every node 
in the GWCC is reachable from every other node in the 
GWCC. As shown in Fig.Q] the GWCC contains the gi- 
ant strongly connected component GSCC, consisting of all 
nodes reachable from each other following directed links. 
All the nodes reachable from the GSCC are referred to 
as the giant OUT component, and the nodes from which 
the GSCC is reachable are referred to as the giant IN 
component. The GSCC is the intersection of the giant 
IN- and OUT-components. All nodes in the GWCC, but 
not in the IN- and OUT-components, are referred to as 
"tendrils" . 



III. COMPONENTS IN TIME- VARYING 
GRAPHS 

In this paper, we consider time-varying graphs, which 
are graphs characterized by links that appear and dis- 
appear over time among a fixed set of nodes. A time- 
varying graph can be described as an ordered sequence 
of graphs, i.e., an ordered set {Gi, G2, ■ ■ ■ , Gm} 01 M 
graphs defined over TV nodes, where each graph G m in 
the sequence represents the state of the network, i.e., the 
configuration of links, at time t m , where m = 1, . . . , M . 



In this notation, the quantity tia ~ ti is the temporal 
length of the observation period. The graphs in the 
sequence can be uniformly distributed over time, i.e., 
t, n +i = t m + At, Vm = 1, . . . , M — 1 HH, or in general 
they can correspond to any ordered sequence of times 
such that ii < £2 < • • • < tu [30l ] . In a more compact 
notation, we denote the graph sequence as Q = G\t u tMY 
Each graph in the sequence can be either undirected or 
directed. Consequently, the time- varying graph Q can be 
described by means of a time-dependent adjacency ma- 
trix A(t m ), m = 1, . . . , M, where aij(t m ) are the entries 
of the adjacency matrix of the graph at time t m . This 
matrix is in general non-symmetric. If we discard the 
time-ordering of the links of a time- varying graph Q and 
consider all links as concurrent in time, we obtain its 
corresponding aggregated static graph. Different ways of 
constructing an aggregated graph from a temporal net- 
work have been indeed proposed. It is possible to assign 
a weight to each edge of the aggregated graph, accord- 
ing to the frequency or the strength of the connection, 
and also to preserve the directionality of edges of the 
original temporal system. However, all these static rep- 
resentations fail to capture the temporal correlation of 
links and, consequently, are not able to represent causal- 
ity relationships. In panel a) of Fig. [5] we report a sim- 
ple time- varying graph ^[f Ij t 4 i with N = 5 nodes and 
M = 4 edges and, in panel b), the corresponding ag- 
gregated static graph (undirected and unweighted) . It is 
worth noticing that the aggregated graph discards most 
of the richness of the original time-varying graph. For 
instance, three paths exist between node 1 and node 5 in 
the static aggregated graph, namely 1-4-5, 1-2-5 and 1-2- 
4-5, while in the time- varying graph there is no temporal 
path from node 1 to node 5. The problem of defining con- 
nectedness and components in time- varying graphs looks 
more similar to the case of directed static graphs than to 
the case of undirected static graphs. In fact, even if each 
graph G m , m = 1, . . . , M in the sequence is undirected, 
the temporal ordering of the graphs naturally introduces 
a directionality. For instance, in the time-varying graph 
G[ti,t 4 ] reported in Fig. [21 there exists a path connecting 
node 5 to node 1 (i.e., the link 052 at time t\ and the 
link 021 at time £3) but there is no path which connects 
node 1 to node 5. An immediate consequence of this fact 
is that node 5 can send a message to node 1 at time t\, 
while node 1 cannot send a message to node 5. 

In order to define node connectedness for a time- 
varying graph, we first need to introduce a mathe- 
matical definition of reachability for an ordered pair 
of nodes i and j. We say that i can reach j if i 
can send a message to j directly or through a time- 
ordered sequence of contacts. In mathematical terms 
this implies the existence of a walk connecting i to j. 
In a time-varying graph, a walk, also called temporal 
walk, from node i to node j is defined as a sequence 
of L edges [(n ro ,n ri ), (n ri ,n r2 ), . . . , (n rL _ 1 ,n rL )}, with 
n ro = i, n rL = j, and an increasing sequence of times 
t ri < t r . 2 < . . . < t rL such that a nr >rh . (r; ) ^ I = 
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(a) 




(b) 

FIG. 2. A time- varying graph Q consisting of a sequence of 
M = 4 graphs with N — 5 nodes (panel a) and its corre- 
sponding aggregated static graph (panel b). The static rep- 
resentation of graphs discards time ordering of links and time 
correlations of paths. In the aggregated graph node 1 and 
node 2 are neighbors, but in the original time- varying graph 
they are directly connected only in one of the four graphs of 
the sequence, namely in Gt 3 ■ Moreover, in the aggregated 
graph a path exists from node 1 to node 5 and vice-versa, 
while in the time-varying graph there exists a temporal path 
from 5 to 1 but there are no temporal paths from 1 to 5. 

1 , . . . , L [3(1 . A path (also called temporal path) of 
a time-varying graph is a walk for which each node is 
visited at most once. For instance, in the time-varying 
graph of Fig. [2 the sequence of edges [(5, 2), (2, 1)] to- 
gether with the sequence of times t\,t$ is a temporal 
path of the graph. This path starts at node 5 at time t\ 
and arrives at node 1 at time £3. Given the definitions of 
temporal walk and path, we can introduce the concepts 
of temporal connectedness (in a weak and in a strong 
sense) for a pair of nodes. 

A node i of a time- varying graph Q\t u tM] ^ s te m P ora ^y 
connected to a node j if there exists in [t\, t^i] a temporal 
path going from i to j. This relation is not symmetric: if 
node i is temporally connected to node j, in general node 
j can be either temporally connected or disconnected to 
i. In the graph S[t 1 ,t 4 ] °f Fig- HI node 5 is temporally con- 
nected to 1 but node 1 is not connected to node 5. For 
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this reason, we introduce the definition of strong connect- 
edness, which enforces symmetry: 

Definition 1 (Strong connectedness) Two nodes i 
and j of a time-varying graph are strongly connected if 

1 is temporally connected to j and also j is temporally 
connected to i. 

Strong connectedness is a reflexive and symmetric re- 
lation, so that if i is strongly connected to j, then j 
is strongly connected to i. However, this definition of 
strong connectedness lacks transitivity, and, therefore, it 
is not an equivalence relation. In fact, if i and j are 
strongly connected and j and / arc strongly connected, 
nothing can be said, in general, about the connectedness 
of i and I. In the example shown in Fig. [2J node 5 and 

2 are strongly connected and also 2 and 1 arc strongly 
connected, but nodes 5 and 1 are not strongly connected, 
since there exists no temporal path which connects node 1 
to node 5. It is also possible to introduce the concept of 
weak connectedness for a pair of nodes. Similarly to the 
case of static directed graphs, given a time- varying graph 
Q, we construct the underlying undirected time- varying 
graph Q u , which is obtained from Q by discarding the 
directionality of the links of all the graphs {G m }, while 
retaining their time ordering. 

Definition 2 (Weak connectedness) Two nodes i 
and j of a time-varying graph are weakly connected if 
i is temporally connected to j and also j is temporally 
connected to i in the underlying undirected time-varying 
graph Q u . 

Also weak connectedness is a reflexive and symmetric 
relation, but it is not transitive. This definition of weak 
connectedness is quite similar, but not identical, to that 
given for directed static graphs. In fact, two nodes in 
Q can be weakly connected even if there is no temporal 
directed path which connects them, but the temporal or- 
dering of links breaks the transitivity so that if i and j are 
weakly connected and j and I are weakly connected, then 
nothing can be said about the weak connectedness of i 
and I. Notice that the definitions of strong and weak con- 
nectedness given above for time-varying graph are con- 
sistent with those given for static graphs, so that if two 
nodes are strongly (weakly) connected in a time- varying 
graph, then they are also strongly (weakly) connected 
in the corresponding aggregated static graph. The vice- 
versa is trivially not true. We are now ready to give 
the definitions of components associated to a node of a 
time- varying graph Q: 

1. The temporal out-component of node i, denoted as 
OUTx(«): is the set of vertices which can be reached 
from i in the time- varying graph Q . 

2. The temporal in- component of a node i, denoted as 
INt(«) 5 is the sc t °f vertices from which i can be 
reached in the time- varying graph Q. 



3. The temporal weakly connected component of a node 
i, denoted as WCCr(*): is t ne set of vertices which 
i can reach, and from which i can be reached, in 
the underlying undirected time- varying graph Q u . 

4. The temporal strongly connected component of a 
node i, denoted as SCCt(«)j is the set of vertices 
from which vertex i can be reached, and which can 
be reached from i, in the time-varying graph Q. 

Differently from the case of directed static graphs, it 
is not possible to define the strongly (weakly) connected 
components of a time-varying graph starting from the 
definition of connectedness for pairs of nodes. As we ex- 
plained above, this is because the relation of strongly 
(weakly) connectedness for couples of nodes is not an 
equivalence relation. For this reason, we give the fol- 
lowing definition of strongly connected component of a 
time- varying graph: 

Definition 3 (Strongly connected component) 

A set of nodes of a time-varying graph Q is a temporal 
strongly connected component of Q if each node of the 
set is strongly connected to all the other nodes in the set. 

Similarly, a set of nodes is a weakly connected compo- 
nent if each node in the set is weakly connected to all 
the other nodes in the set. The definitions of strongly 
and weakly connected components enforce transitivity, 
but the check of strong (weak) connectedness has to be 
directly performed for every couple of nodes. Let us sup- 
pose for instance that we want to verify if the five nodes 
in the graph Q shown in Fig. [2] form a strongly connected 
component. In the static aggregated graph this check 
has O(K) computational complexity, where K is the to- 
tal number of links in the graph. In fact, we have only 
to check that 2, 3, 4 and 5 are connected to I, which 
can be done by a depth first visit of the graph started at 
node 1, since node connectedness is an equivalence rela- 
tion for static graphs, and a component of a node is also 
a component for the whole graph. On the contrary, for a 
time- varying graph we should check the connectedness of 
all the possible couples of nodes, so that a procedure to 
verify that a set of N nodes forms a strongly connected 
component has computational complexity 0(N 2 ) for ev- 
ery check, instead of 0(K). Such a difference in compu- 
tational complexity is relevant for graphs obtained from 
real world systems, which are usually sparse. Moreover, 
while static directed graphs admit only one partition into 
strongly connected components, for a time- varying graph 
there exists in general more than one possible partition, 
as we shall see in the next section. 



IV. THE AFFINE GRAPH OF A 
TIME- VARYING GRAPH 

In this section we show that the problem of finding 
the strongly connected components of a time-varying 
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graph is equivalent to the well-known problem of finding 
the maximal-cliques of an opportunely constructed static 
graph [HI. We call such a static graph the affine graph 
corresponding to the time- varying graph. It is defined as 
follows: 

Definition 4 (Affine graph of Q) 

Given a time-varying graph Q = Q[ti,t M ]i ^ e asso- 
ciated affine graph Gg is an undirected static graph with 
the same nodes as Q , and such that two nodes i and j 
are linked in Gg if i and j are strongly connected in Q. 

In practice, the affine graph of a time- varying graph can 
be obtained by computing the temporal shortest paths 
between any two pairs of nodes, and then adding a link 
between two nodes i and j of the affine graph only if 
the temporal distance from i to j and the temporal dis- 
tance from j to i are both finite. Another method to 
construct the affine graph is based on the usage of the 
out-components of all the nodes. We start by consider- 
ing the out-component of the first node (let us say i = 1) 
and then we check if for each node j E OUTt{i),J > i 
then also i E OUTt{])- If this is true, we add a link 
between i and j in the affine graph. We then repeat this 
procedure for the second node, i = 2, for the third node, 
i = 3 and so on. We obtain the affine graph by iterat- 
ing over the out-components of all the nodes. It is worth 
noticing that the definition of affine graph is independent 
from the particular definitions of reachability and tempo- 
ral path given here, since it is based only on the existence 
of mutual temporal paths between pairs of nodes. There- 
fore, given any of the different definitions of time- varying 
graph and the corresponding notions of walk and path, it 
is always possible to construct the affine graph associated 
to that particular representation of the system. 

We also observe that the concept of affine graph is 
somehow related to existing concepts of path graphs, 
reachability graphs and influence digraphs pBl [29J. In 
fact, in a path graph (or in an influence digraph) a di- 
rected link exists from node i to node j if there is a tempo- 
ral path from itoj, i.e., if it is possible to reach j starting 
from i. Therefore, we can construct an affine graph from 
the path graph or from the influence digraph associated 
to a time- varying graph, by replacing each couple of re- 
ciprocated directed links with a single undirected edge, 
and removing all the remaining links. 

In Fig. [3] we report the affine graph corresponding to 
the time varying graph shown in Fig. [2j In this graph, 
node 1 is directly connected to nodes {2,3,4}, since it 
is temporally strongly connected to them in the time- 
varying graph. Similarly, node 2 is connected to nodes 
{1,3,4,5}, node 3 is connected to {1,2}, node 4 is con- 
nected to {1,2,5} and node 5 is connected to {2,4}. 
Hence, the affine graph Gg has only 7 of the 10 pos- 
sible links, each link representing strong connectedness 
between two nodes. 

We briefly report here some definitions about graph 
cliques. Given an undirected static graph, a clique is a 
complete subgraph, i.e., a subgraph in which all the nodes 




FIG. 3. The affine graph Gg associated to the time-varying 
graph Q reported in Fig. [2] The affine graph is static and 
undirected, and each of its maximal-cliques corresponds to 
a strongly connected component of the original time-varying 
graph g. 

are directly linked to each other. A maximal- clique is a 
clique that is not included in any larger clique, while a 
maximum- clique is a maximal- clique whose size is equal 
to or larger than those of all the other cliques [39| . 

By construction, a clique of the affine graph Gg, con- 
tains nodes which are strongly connected to each other, 
so that the maximal- cliques of the affine graph, i.e., all 
the cliques which are not contained in any other clique, 
are temporal strongly connected components (SCCt) 
of Q. Similarly, all the maximum- cliques of the affine 
graph Gg, i.e., its largest maximal-cliques, are the largest 
temporal strongly connected components (LSCCt) of Q. 
Therefore, the affine graph can be used to study the 
connectedness of a time-varying graph, and the prop- 
erties of the strongly connected components of a time- 
varying graphs can be obtained from known results about 
maximal-cliques on static graphs. For instance, the prob- 
lem of finding a partition of Q that contains the mini- 
mum number of disjoint strongly connected components 
is equivalent to the well-known problem of finding a par- 
tition of the corresponding affine graph Gg in the small- 
est number of disjoint maximal-cliques [38(. Unfortu- 
nately, this problem is known to be NP-complete, and in 
practice can be exactly solved only for small graphs. In 
the case of the affine graph in Fig. [31 it is possible to check 
by hand that there are only three possible partitions of 
Gg into maximal-cliques, namely: 

1. {1,2,3}(J{4,5} 

2. {1,2,4}U{3}U{5} 

3. {2,4,5}U{1,3} 

Notice that the second partition contains two isolated 
nodes, which are indeed degenerated maximal-cliques. 
Therefore, the original time-varying graph admits only 
two different partitions into a minimal number of non- 
degenerated strongly connected components, namely into 
two components containing at least two nodes each. One 
possible partition of our network Qu >t4 ] is composed of 
the components {1, 2, 3} and {4, 5}, while the other par- 
tition consists of {2,4,5} and {1,3}. If we discard the 
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temporal ordering of links, we obtain different results. 
In fact, the aggregated static graph shown in Fig. [2] has 
only one connected component, which includes all the 
five nodes. 

Other interesting results stem from the mapping into 
affine graphs and from the following well known results 
for cliques in graphs. 

1 . Checking if a graph contains a clique of a given size 
k has polynomial computational complexity, and 
precisely 0(N k k 2 ) go[. 

2. The clique decision problem, i.e., the problem of 
testing whether a graph contains a clique larger 
than a given size k, is NP-complete [38|. There- 
fore, any algorithm that verifies if a time-varying 
graph has a strongly connected component whose 
size is larger than a fixed value k, has exponential 
computational complexity 

3. Listing all the maximal-cliques of a graph has expo- 
nential computational complexity, namely 0(3 N / 3 ) 
on a graph with N nodes [IH [42[ . Consequently, 
finding all strongly connected components of a 
time- varying graph with A" nodes, requires an 
amount of time which exponentially grows with N. 

4. The problem of finding a maximum-clique for 
an undirecte d grap h is known to be hard-to- 
approximate [43l445| . and an algorithm that finds 
maximum-cliques requires exponential time. This 
means that, if P ^ NP, there exists no polyno- 
mial algorithm to find approximated solutions for 
the maximum-clique problem. The best known al- 
gorithm works in 0(~ 1.2^) for a graph with N 
nodes 0,113. 

5. The problem of determining if a graph can be parti- 
tioned into K different cliques is NP-complete, and 
consequently also the problem of finding the mini- 
mum number of cliques that cover a graph, known 
as the minimum clique cover, is NP-complete [38j . 
This means that there exists no efficient algorithm 
to find a partition of a time-varying graph made 
by a set of disjoint strongly connected components. 
Moreover, there are in general more than one parti- 
tion of a graph into maximal-cliques, so that a time- 
varying graph cannot be uniquely partitioned into 
a set of disjoint strongly connected components. 

The existence of a relation between the strongly con- 
nected components of a time-varying graph and the 
maximal-cliques of its affine graph implies that it is prac- 
tically unfeasible to find all the strongly connected com- 
ponents of large time-varying graphs. The problem can 
be exactly solved only for relatively small networks, for 
which it is computationally feasible to enumerate all the 
maximal-cliques of the corresponding affine graphs. Even 
if, in many practical cases, it is possible to find only the 
maximal-cliques up to a certain size k, we can still obtain 



some information about the maximum value of k to be 
checked. First of all, in order to have a clique of size k 
the graph should have at least k nodes having at least k 
links. Moreover, each clique of order k > 3 has exactly 
(g) sub-cliques of order 3, so that in order for a subgraph 
to be a clique of order k, the graph should have at least 

(3) triangles. This means that there is a relation between 
the number of triangles of the affine graph and the size 
of its maximum-cliques. In particular, the number of ex- 
isting triangles in the affine graph sets an upper bound 
for the size of the largest admissible maximal-cliques of 
the graph. 



V. RESULTS 

In this section we extract and analyze node and graph 
components of time-varying graphs constructed from 
three different data sets of temporal social networks. The 
first is the data set of human interactions produced by 
the Reality Mining Project [48[ , which was a large exper- 
iment devised and performed by researchers at the Me- 
diaLab, Massachusetts Institute of Technology, US. One 
of the aims of this experiment was to record the contacts 
among students, staff and faculty members at the Mas- 
sachusetts Institute of Technology The resulting data 
set contains co-location information among 100 individ- 
uals during six months, from the end of June 2004 to the 
end of December 2004, sampled by means of Bluetooth- 
enabled mobile phones at regular intervals of five min- 
utes. The second data set contains all the contacts among 
78 people attending the INFOCOM 2006 conference 
Similarly to the Reality Mining data set, also the INFO- 
COM 2006 dataset has been constructed using records of 
scans made by Bluetooth devices provided to conference 
participants, but covers a shorter period of time (only 
five days). The third data set contains the records of 
communications among ~ 100.000 Facebook users living 
in Santa Barbara (CA), occurred over a period of six 
months, from January 2009 to June 2009 [5(J. In the 
following we denote the Reality Mining, INFOCOM and 
Facebook data sets as RM, IC and FB respectively 

When contacts are inferred from Bluetooth scans, we 
make the reasonable assumption that two individuals are 
co-located, i.e., they are at the same place, at a given 
time, if their respective devices detect each other pres- 
ence. In fact, Bluetooth devices are able to detect simi- 
lar devices within a very limited range, usually between 5 
and 10 meters. At each time t, a co-location graph can be 
obtained by connecting through undirected links all the 
nodes which are co-located at that time. For the RM 
and the IC data sets we constructed several time- varying 
graphs, made of sequences of co-location graphs obtained 
at regular intervals of 5 minutes. For the dataset of FB 
communication, the concept of co-location is not mean- 
ingful, since the communication is usually asynchronous 
and mediated by the Internet. In this case, we construct 
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a contact graph for each hour, where a link between two 
nodes does exist if the nodes have exchanged at least one 
message in the corresponding interval. We have chosen 
these three data, sets because they represent different typ- 
ical aspects of temporal social networks. First of all, RM 
provides a good example of recurrent face-to-face inter- 
actions, like those we experience in our everyday life. At 
the same time, IC is a relevant instance of bursty face- 
to-face interactions, which are typical in gatherings and 
meetings. Finally, FB is the prototypal model of on-line 
communication, which is different from face-to-face inter- 
action because it does not require physical co-location to 
be facilitated. 

We start analyzing the distribution of the sizes of the 
node temporal in- and out- components in RM, by con- 
sidering time- varying graphs obtained from the first half 
and the second half of a week. In particular, in Fig. [4] we 
consider week 11, which roughly corresponds to the be- 
ginning of the Fall term. For each node, we report the size 
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FIG. 4. Size of the temporal in-component (panel a) and 
out-component (panel b) for each of the N — 100 individuals 
during week 11 of the RM data set. Red circles and blue 
squares correspond, respectively, to the beginning of the week 
(WB) and to end of the week (WE). For comparison, the 
size of the largest connected component of the corresponding 
aggregated static graph are reported as dashed red line (WB) 
and solid blue line (WE), respectively. 

of its temporal in- and out-component during the begin- 
ning of the week (WB), namely from Monday 12:00am 
to Thursday 11:59am (red circles), and during the end 
of the week (WE), namely from Thursday 12:00pm to 
Sunday 11:59pm (blue squares). As shown in the figure, 
during WB almost all the nodes have temporal in- and 
out-component of size 72. Conversely, during WE, we 
observe a wider distribution of the sizes of temporal in- 
and out-components. In particular, in panel (a) we no- 
tice a group of nodes having an in-component of size 53, 
another group whose in-component contains around 40 
nodes, and other nodes with in-component of size smaller 



than 30. Similarly, in panel (b), there is a group of nodes 
whose out-component contains around 60 nodes, a sec- 
ond group of nodes with out-component sizes between 
40 and 50, and many other nodes having out-component 
with less than 40 nodes. The observed small variability 
in the size of node components during WB, is due to the 
fact that students and faculty members have more op- 
portunities to meet and interact at lectures during WB. 
Even if not all students attend the same classes, and not 
all professors teach to all the students, there is a high 
probability that two individuals would be connected by 
longer temporal paths. Conversely, during WE, the stu- 
dents usually meet other students in small groups, and 
they usually do not meet professors and lecturers, except 
for the classes held on Thursday afternoon and on Friday. 
As a result, the size of the in- and out-components during 
WE exhibits large fluctuations from node to node. Such 
fluctuations are lost in a static graph description, which 
aggregates all the links independently of their time or- 
dering. In fact, the static aggregated graphs correspond- 
ing respectively to WB and WE, have only one giant 
connected component which contain the majority of the 
nodes, while the remaining nodes are isolated. As com- 
parison, the size of the giant component of the aggre- 
gated static graphs for WB and WE are also reported in 
Fig.|H respectively as dashed red line and solid blue line. 
Notice that the static aggregated graph corresponding 
to a co-location time- varying graph is intrinsically undi- 
rected. Therefore, the in- and out-components of a node 
in this graph coincide and correspond to the component 
to which the node belongs. Moreover, in a static ag- 
gregated graph all the links (and consequently also all 
the paths) are always available, so that all the nodes in 
the same connected component have the same compo- 
nent size. As a result, the variability in the node con- 
nectedness of the time- varying network, which is evident 
from the distribution of circles and squares in Fig. 21 
is flattened down in the aggregated static graphs. In 
the latter case, all information about network connect- 
edness is represented by a single value, namely the size 
of the largest connected component, which does not pro- 
vide any information about the mutual reachability of 
two generic nodes of such a component. In particular, 
the size of the giant connected component of the static 
aggregated graph is equal to 74 during WB and to 66 
during WE, despite the fact that in the latter interval 
the majority of nodes has much smaller temporal in- 
and out-components. Similar results are obtained for the 
temporal in- and out-components of nodes in IC and FB. 
However, we do not report here the corresponding figures 
for brevity. 

We now analyze the strongly connected components 
of different time- varying graphs obtained from the three 
data sets under consideration. In particular, we focus on 
the Fall term of RM (namely weeks from 10 to 19 in the 
original data set), the last four days of IC (namely days 
2 to 5), and the first twelve weeks of FB (from January 
to the end of March 2009). In Table Q] we report the rel- 
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evant structural properties of different affinc graphs con- 
structed from RM. In particular, we consider the time- 
varying graphs constructed in the first 24 hours (Mon- 
day) of ten consecutive weeks (from week 10 to week 19), 
and the time-varying graphs corresponding to the whole 
weeks. At the scale of one day, we observe large fluctu- 
ations in the measured values. The number of links K 
ranges from 105 in week 12 to 1485 in week 15, while the 
number of triangles T is in the range [307, 22096], with a 
mean value around 10000 and a standard deviation equal 
to 6932. This variance is due to the fact that, even if the 
daily activity of each individual is, on average, almost 
periodic, in a particular day we can observe a peculiar 
temporal pattern of connections. For instance, some stu- 
dents decide to skip a class or the lessons are suspended 
for public holidays. This is exactly what happens on 
week 12. Monday of week 12 is September II th 2004, and 
corresponds to the Patriot Day, a national holiday intro- 
duced in the US in October 2001, designated in memory 
of the 2977 killed in the September 11 th , 2001 attacks. 
Therefore, we observe the minimum connectivity and the 
minimum number of triangles on week 12 of RM, because 
all teaching activities were suspended, and students did 
not participate to lessons as usual. At the scale of one 
week, instead, the number of links and the number of tri- 
angles T are much more stable over time: K has a mean 
value of 2717 and a standard deviation of 342, while T 
has a mean value around 62000 and a standard devia- 
tion around 12000. For the graphs of Mondays, we find 
that the number N s and the average size (s) of maximal 
cliques of the affinc graphs fluctuate over the ten weeks. 
In particular, we observe relative smaller values of N s and 
(s) during Mondays of weeks 10 to 14 than in weeks 15 
to 19, which is probably due to the relatively lower num- 
ber of links and triangles. Conversely, if we consider the 
size S of the largest strongly connected component (i.e., 
the largest maximal-clique of the affine graph), we no- 
tice that it is not strongly correlated with K and T. For 
instance, the size of the largest strongly connected com- 
ponent found at Monday of week 11 (S = 29) is equal to 
that observed at Monday of week 16. However, in week 
11 the affine graph has a much smaller number of links 
and triangles than in week 16. Moreover, on Monday of 
week 14 we have a maximal-clique of size 27, even if the 
number of links and triangles is higher than on Monday 
of week 1 1 . 

At the scale of a week, we still observe relevant fluc- 
tuations of average size (s), in the range [39.3,62.5], but 
the number of strongly connected components is stable 
around N s = 10, with a maximum of Ns = 15 in week 14. 
These results confirm that the size of the largest strongly 
connected component of a time-varying graph is deter- 
mined by the actual configuration of links and triangles 
of the corresponding affine graph, and not only by their 
relative number, and in general depends on the temporal 
scale at which we observe the system. 

Examining the results reported in Table HI we notice 
that the affine graphs for Mondays of all the weeks of 
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TABLE I. Structural properties of the affine graph corre- 
sponding to the time- varying graph of the first 24 hours of the 
week (Monday), and to the whole weeks in the Fall term of 
RM. We report the number of links (K), number of triangles 
(T), number of maximal cliques (N 3 ), average size of maxi- 
mal cliques ((s)), size of the largest maximal clique (5*), num- 
ber of largest maximal cliques (Ns), number of nodes in the 
union (Nu) and in the intersection (Ni) of all largest maximal 
cliques. The size of the giant component of the correspond- 
ing static aggregated graph (C) is reported in the rightmost 
column. 







K 


T 


N s 


is) 


S 


N S 


Nu 


Ni 


C 




6:00 


376 


2061 


17 


8.1 


19 


2 


20 


18 


44 




7:00 


793 


7514 


16 


15.8 


29 


1 


29 


29 


50 




8:00 


1922 


35472 


13 


35.8 


56 


2 


57 


55 


69 




9:00 


1825 


35598 


4 


29.5 


60 


1 


60 


60 


61 




10:00 


2252 


48459 


6 


50.8 


64 


1 


64 


64 


68 




11:00 


2142 


45510 


3 


43.0 


65 


1 


65 


65 


66 


>-< 




12:00 


2346 


52394 


1 


69.0 


09 


1 


69 


69 


69 


o 
— 


13:00 


2621 


61270 


2 


69.0 


72 


1 


72 


72 


73 


O 


14:00 


2389 


53297 


2 


56.5 


69 


1 


69 


69 


70 


15:00 


2330 


51382 


4 


63.5 


66 


1 


66 


66 


69 




16:00 


1991 


40153 


5 


46.6 


58 


3 


60 


56 


65 




17:00 


1862 


36165 


5 


44.2 


58 


1 


58 


58 


62 




18:00 


1595 


27979 


6 


39.3 


51 


1 


51 


51 


58 




19:00 


1177 


18424 


2 


25.5 


49 


1 


49 


49 


49 




20:00 


1110 


16393 


8 


16.0 


46 


1 


46 


46 


49 






K 


T 


N s 


<*> 


.5 


N s 


Nu 


N r 


C 


CO 


2 


2920 


72638 


4 


55.8 


75 


1 


75 


75 


78 


if 


3 


2776 


67435 


3 


50.3 


74 


1 


74 


74 


76 




4 


2225 


47287 


5 


44.0 


63 


1 


63 


63 


68 


u 


5 


1602 


27617 


9 


31.7 


53 


1 


53 


53 


60 



TABLE II. Structural properties of the affine graphs corre- 
sponding to time- varying graphs of different hours of the third 
day and of each of the four days of IC. The graph correspond- 
ing to each hour includes all the contacts recorded in that 
hour, so that, for instance, the first graph is constructed from 
the interactions observed from 6 : 00 to 6 : 59. Legend as in 
Tabled] 

RM admit a single LSCCt, except in week 19 where two 
LSCCys of size S — 38 emerge. For this reason, we also 
looked at the number of nodes Njj which participate to at 
least one LSCCt, and at the number Nj of nodes which 
participate to all LSCCts. These numbers correspond, 
respectively, to the number of nodes found in the union 
and in the intersection of all LSCCts. An interesting 
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result is that Nj = 34 on Monday of week 19, so that 
34 nodes participate to both maximal 42-node cliques. 
These 34 nodes play a very important role in the struc- 
ture of the network. If we remove just one of them, then 
the resulting affinc graph does not have a clique of size 
42 any more, and consequently the size of the LSCCt of 
the remaining time- varying graph is smaller than 42. At 
the same time, removing all these Ni nodes will cause a 
significant reduction in the size of LSCCts, in the num- 
ber of triangles of the affine graph and, consequently, in 
the number of SCCtS. The nodes that participate to at 
least one LSCCt are important for the diffusion of in- 
formation throughout time-varying graphs. In fact, if a 
message is passed to one of these Nu nodes, then there 
is a high probability that the message could be delivered 
to all the other Njj — 1 nodes in the union of all LSCCt- 
Similar results arc observed also at the scale of a week, 
where almost all graphs have just one LSCCt, except for 
week 18 and 19 where we found N$ = 2 and Ns = 3, 
respectively. 

In the rightmost column of Table U we report, for each 
time- varying graph, the size C of the giant component of 
the corresponding static aggregated graph. At the scale 
of one day, the value of C is always much larger than S, as 
a consequence of the fact that the static representation 
of the time-varying graph systematically overestimates 
node connectedness and paths availability. Conversely, 
at the scale of one week the size of the LSCCt is com- 
parable with the size of the largest strongly connected 
component of the aggregated graph. In fact, when we 
consider longer time periods, longer temporal paths ap- 
pear and, consequently, there is a higher probability that 
a large number of pairs of nodes become mutually reach- 
able. The affine graphs of entire weeks of RM are almost 
complete graphs and look similar to each other, so that 
we observe a relatively high correlation between S and C 
at the scale of one week, even if there is no correlation be- 
tween S and C at the scale of a single day. For instance, 
on Monday of week 16 we observe the maximum value of 
C, namely C = 70, while the corresponding time- varying 
graph has a largest strongly connected component of size 
S = 29, which is relatively small compared to the other 
weeks. Conversely, on Monday of week 13 we observe 
a relatively small giant component, with C — 59 nodes, 
while the size of the largest strongly connected compo- 
nent at the same day is S = 36. 

In Table [TT] we consider the time- varying graphs con- 
structed at different hours of the day 3 of IC, and those 
obtained for each of the four days of the conference. In 
this case, the system at the scale of one day is satu- 
rated: there arc just a few strongly connected compo- 
nents and only one LSCCt, which roughly correspond 
with the largest connected component of the aggregated 
graph. In fact, the dynamics of interactions between par- 
ticipants at a conference differs a lot from the recurrent 
daily interactions which arc typical in RM. Usually, dur- 
ing a single day of conference the participants attend 
different plenary sessions, and consequently they have 



many opportunities to be co-located with a large num- 
ber of other participants for a relatively long period (at 
least a couple of hours). Therefore, all the participants 
to the same session usually form a single strongly con- 
nected component. However, if we observe the system at 
the scale of one hour, we notice that there are relevant 
fluctuations in the number of links and triangles, and in 
the number and size of strongly connected components. 
In the graphs constructed from contacts recorded before 
8:00 we observe a relatively smaller number of links and 
triangles and smaller values of S. The number of links 
remains almost stable until 17:00, while S steadily in- 
creases up to S = 72 at 13:00 (lunch time) when the 
size of the LSCCt reaches a maximum. Then, both the 
number of triangles and the size of the LSCCt decrease, 
and the graph at 20:00 has S = 46 and T = 16393. Also 
at the scale of one hour we observe a clear correlation 
between S and C, Ns = 1 and the average size (s) of 
SCCtS is comparable with the size of the largest strongly 
connected component. 
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119 


119 


44464 


6 


345943 


4805124 


392966 


31.0 


91 


1 


91 


91 


44117 


-a 


7 


346856 


5527618 


149156 


23.0 


145 


1 


145 


145 


44107 




8 


313564 


3176483 


195238 


25.2 


138 


1 


138 


138 


44886 


FB 


9 


493502 


14700802 


1160920 


66.7 


149 


1 


149 


149 


46790 


10 


282759 


2966668 


154016 


21.5 


87 


1 


87 


87 


43693 




11 


217796 


2277735 


29077 


13.2 


103 


2 


104 


102 


39950 



TABLE III. Structural properties of the affine graphs corre- 
sponding to the time-varying graphs of weeks (upper rows) 
and pairs of adjacent weeks (lower rows) of FB data set. Leg- 
end as in Table [I] 

Finally, in Table lnTl we report the structural properties 
of the affine graphs constructed from FB. We considered 
12 time- varying graphs at the scale of a week and 11 
time- varying graphs at the scale of 14 days, i.e. couples 
of subsequent weeks. In this case the number of nodes 
of the affine graphs was quite large, so we used HPC fa- 
cilities provided by the Edinburgh Parallel Computing 
Center (EPCC), in order to run a parallelized version of 
the algorithm described in [42| on 128 processors. At the 
scale of a week, we observe significant fluctuations in the 
structural properties of the affine graphs corresponding 
to distinct weeks: there is a high variability in the num- 
ber of links (K is in the range [43400, 79645]) and in the 
number of triangles (T is in the range [32655,200061]). 
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However, the size of the LSCCt is relatively small, com- 
pared to the number of nodes in the graph. In fact, S 
remains in the range [16,41], which is quite surprising 
for a graph with more than 100000 nodes. At the same 
time, there is no correlation at all neither between S and 
C nor between S and T. A similar result is obtained 
when considering pairs of adjacent weeks: in this case, 
the size of the LSCCt is slightly larger (the maximum is 
observed for the 9 th pair of weeks, when S = 149), but 
S is still two orders of magnitude smaller than C . The 
case of FB is emblematic of the implicit weaknesses of a 
component analysis based on static aggregated graphs. 
According to traditional static graphs theory, the giant 
component of a network constructed from two adjacent 
weeks of Facebook communication contains at least 32% 
of the 100000 nodes, while a temporal analysis reveals 
that less than 0.15% of the nodes are indeed mutually 
reachable. 



VI. CONCLUSIONS 

Conventional definitions of connectedness and compo- 
nents proposed so far have only considered aggregate, 
static topologies, neglecting important temporal infor- 
mation such as time order, duration and frequency of 
links. In this work we have extended the concepts of 
connectedness to the case of time-varying graphs, and 
we have introduced definitions of node and graph compo- 
nents which take into account duration, times of appear- 
ance and temporal correlations of links. The proposed 



temporal measures are able to capture variations and 
fluctuations in the linking patterns, typical of many real 
social and biological systems. As a first application, we 
have studied three databases of human contacts, show- 
ing that variations in the pattern of connections among 
nodes produce relevant differences in the size and num- 
ber of temporal strongly connected components. We have 
pointed out the important role played by nodes that be- 
long to different strongly connected components at the 
same time, and we have also analyzed how the size of 
temporal strongly connected components depend on the 
observation interval. We hope that our formalism will 
be useful to analyze other data sets of time- varying net- 
works that will be available in the near future, and to 
better characterize dynamical processes that take place 
on these networks, such as diffusion of information and 
spreading of diseases. 
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